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ABSTRACT 

We have constructed one-phase quasi-periodic solutions of the curve equation induced 
by the mKdV equation. The solution is expressed in terms of the elliptic functions of 
Weierstrass. This solution can describe curve dynamics such as a vortex filament with axial 
velocity embedded in an incompressible inviscid fluid. There exist two types of curves (type- 
A, type-B) according to the form of the main spectra of the finite-band integrated solution. 
Our solution includes various filament shapes such as the Kelvin-type wave, the rigid vortex, 
plane curves, closed curves, and the Hasimoto one-solitonic filament. 
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1 Introduction 



Space curves are of interest in mathematics and physics. This is especially true with how 
these curves' mathematical structures are related to integrable equations. The first ap- 
pearance of such a relationship was between the sine-Gordon equation and the differential 
geometry for surfaces of constant negative curvature [1]. Another interesting development 
was from Hasimito who showed the motion of vortex filament motion is related to the nonlin- 
ear Schrodinger equation [2]. Motivated by this development, Lamb studied various helical 
curves that were described by well-known integrable equations [3]. He showed that there 
exists a curve that has a constant torsion and its curvature is described by the modified 
Korteweg-de Vries (mKdV) equation. 

A special class zero torsion curves are the planar curves, which already appeared in the 
famous old problem of elastica. The elastica problem was solved by Euler. Typical solutions 
of the elastica problem, found by a numerical method, can be seen in the book [4], which 
are essentially those solutions of Euler. After that, this problem has occurred repeatedly 
in diverse forms. For example, Ishimori [5] found that the soliton solution of the mKdV 
equation is essentially the same as the loop soliton solution found in [6]. The loop soliton 
described the loop solitary wave propagating along a stretched rope. Mumford showed that 
computer vision problems were related to the mKdV equation [7] . A new form of the solution 
was given by Mumford which was expressed by the theta functions of the genus-one Riemann 
surface [7]. More general solutions, corresponding to the genus-N Riemann surfaces, were 
given by Matsutani in terms of the Weierstrass sigma functions and their generalizations [8] . 
On the other hand, [9] and [10] related the modified Korteweg-de Vries (mKdV) equation to 
motions of curves in a plane. 

The above development of the elastica problem was confined to space curves having zero 
torsion and laying on a plane. On the other hand, the problem posed by Lamb, finding 
space curves related to the mKdV equation having a constant torsion and thus moving in 
three-dimensional space, was not investigated thoroughly in the research literature. As we 
can see in Section 2.4, this problem is related to the motion of a vortex filament with axial 
velocity and moving in an incompressible inviscid fluid [11]. Generally, constructing a curve 
having a given curvature and torsion is a difficult problem. It can be carried out for a simple 
from of curvature and torsion only [12]. For example, curves having a curvature of the 
one-soliton and constant non-zero torsion are explicitly constructed in the Appendix of [3]. 
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Curves, having a more complex form of curvature, have been constructed for special cases 
only, by applying methods specially adopted to the given problem [13]. 

In our work, described in this paper, we construct space curves whose curvature is de- 
scribed by the quasi-periodic solution of the mKdV equation and curves which have constant 
non-zero torsion. The constructed curve is of the Kida-type stationary solution [14]. Equa- 
tions that describe this Kida-type stationary solution were already known in [11], but were 
not solved because of their complicated form. The present calculation applied the method 
of the "modified finite-band integration". The modified method was invented to solve the 
so-called effectivization problem, which is related to the extraction of the real solutions from 
the finite-band integration method [15, 16]. This modified method gives solutions in the 
simplest but important one-phase case. The modified method gave physical solutions in a 
rather simple form , which made this method useful to apply to real situations. Recently it 
was used to obtain quasi-periodic solutions of various integrable theories expressed by the 
Weierstrass elliptic functions [17, 18]. It was then applied to space curve problems, such as 
the Heisenberg model, filament motions under the localized induction approximation, and 
the Lund-Regge vortex motions [19, 20, 21]. These problems were related to the nonlinear 
Schrodinger equation and complex sine-Gordon equation, respectively. 

A recent related calculation was the application of the Sym-Tafel formula to the problem 
of filament motions under the localized induction approximation [22]. Ref. [22] also ob- 
tained real quasi-periodic solutions expressed in terms of the Riemann theta function. This 
approach, based on the standard finite-band integration method, has the merit of obtaining 
more general solutions of multi-phase case. But these solutions through the higher-genus 
Riemann surfaces have complicated form, including numerical integrations to explictly enu- 
merate parmeters appearing in the solutions. Thus it was difficult for physicists to apply 
the quasi-periodic solutions of these type to the analyses of experimental results. This was 
the reason of appearing one-phase solutions of various integrable equations in the previous 
paragraph by using the modified finite-band integration methods. Note that there generally 
exist multi-phase solutions of various integrable equations in terms of the Riemann theta 
function, though they are not widely used in the analyses of physical problems [17]. In fact, 
the present paper gives solutions expressed in terms of the Weierstrass' elliptic fuctions with 
explicit parameters. Using the known properties of the Weierstrass' functions, we can trace 
the variation of physical properties, like moving or rotation velocities, of the solutions along 
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with the spectral parameters. 

The derivation of the curve configurations, described in the present paper, is similar to 
the previous cases of [19, 20, 21]. But there exists an important difference with the present 
calculation as compared to the previous ones, related with imposing conditions such as the 
constancy of torsion r. It then results in a constraint on the main spectra and classifies 
the calculated curves into two categories (type- A and type-B). This restriction on the main 
spectra is not observed in the previous calculations of the modified finite-band integration 
method [17, 20, 21]. Solutions of a more simple type, discussed in [11], including the one- 
solitonic curve, the helicoidal vortex filament, and the Euler's elastica, are obtained by the 
reduction from the curve of the present calculation. 

The general setup for the curve problem of the mKdV equation is given in Section 2. 
Explicit construction of the one-phase periodic solution is followed using a modified version 
of the finite-band integration method in Section 3. Especially, imposing special conditions of 
the mKdV problem is given in Section 3.4. The resulting formulae involve Weierstrass' elliptic 
functions. In Section 4, various specific solutions from the reduction of periodic solutions 
are explained using the functional relations of Weierstrass' functions. Some numerical plots, 
as well as an explicit check of solutions, are accomplished with the help of the symbolic 
package, Mathematica. Section 5 follows with a discussion. 

2 mKdV-type space curves 
2.1 Equations of curve motion 

The space curve in three dimension is described by {t, n, b} vector (the unit tangent, normal, 
and binormal vectors, respectively). Here, t = dr/dz, where r(z,z) is the position vector of 
a point on the curve parametrized by the arc length z at time z. They satisfy the following 
Serret-Frenet equation, 



where d = d/dz and k, t are respectively the curvature and torsion. Generally, motion of a 
point on the curve can be specified in the form 



dt = ten, <9n = —Kt + rb, db = — rn, 



(1) 



r = dr = Un + Vb + Wt, 



(2) 
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where d = d/dz is the time derivative. Especially, the mKdV-type curves with a constant 
torsion r is described by [3] 

U = -Ok, V = 2kt, W = 3r 2 - k 2 /2. (3) 

The family of curves with r = are shown to be related with the problem of arc length 
conserving (nonstretching) motion and/or area conserving motion on a plane [9, 10]. It 
was first analyzed by Euler in his study of the thin rod, the so-called the elastica problem. 
Various possible shape of curves are illustrated in [4], including the famous figure-8 shape. 
These type of solutions are obtained by assuming that curves lye on a plane, thus leading to 
curves of zero torsion. On the contrary, the interesting problem of non-zero torsion has not 
been dealt in depth in literature and only the simple solutions has been obtained [3]. 

2.2 Integrable structure of the mKdV equation 

We first start with the zero curvature condition of the mKdV equation, 

[d + U KV , d + V KV } = 0, (4) 

where 

/ -{d 2 n)/2- k 3 /A\ . (-k 2 /A (9k)/2\ 

Vkv ~ {(d 2 K)/2 + K 3 /A ) +lT \(dK)/2 k 2 /A J 

The mKdV variable k = k(z, z) becomes the curvature of the corresponding curve problem, 
while the spectral parameter r becomes the torsion. Note that r is a constant and it does 
not depend on z, z. The zero curvature condition in Eq. (4) is the compatibility condition 
for the overdetermined system of the following associated linear equations 

(d + U KV )$ = 0, {d + VKv)& = 0. (6) 

By explicitly enumerating the zero curvature condition at the level of 0(r°), we obtain the 
mKdV equation, 

3k + d 3 K + ^-k 2 8k = 0, (7) 
while other levels of 0(r % ),i = 1,4 just becomes identities. 
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2.3 The Sym-Tafel formula 

Now we introduce the duality relation between the curve problem and the mKdV system. 
We write down the curve position vector r = (ri,r 2 ,r 3 ) in terms of $ in Eq. (6) and shows 
that the linear equation in Eq. (6) implies the Serret-Frenet equation in Eq. (1). For this, 
we introduce the following Sym-Tafel formula, 

r = 2i$" 1 T f-$, (8) 

OT 

where r = ^ r i°"j^ — 1 ; 3, with the Pauli matrices Oi [23, 24]. (From now on, we will omit 
the summation notation.) 

To check that Eq. (8) gives the Serret-Frenet equation, we first notice that 

(9 f) BU 

dr = -2i$- 1 a$$" 1 — $ + 2i$" 1 — (0$) = -2i^— — $ = $~V 3 $. (9) 

OT OT OT 

In other words, the unit tangent vector t = dr = (ti, t 2 , t 3 ) is given by 

t = ti<Ti = dr t a t = dr = $-V 3 $. (10) 

Note that UU = |Trt 2 = I. This equation shows that the tangent vector t is given by 
the rotation of k (the unit vector along the z-axis), where the rotation is induced by the 
similarity transformation of $. In a similar way, we find that 

3t = dtiVi = d 2 r = -[^d®, $~V 3 $] = ^[Ukv, <t 3 ]$ = -k®' 1 ^®. (11) 

Eq. (11) is one of the Serret-Frenet equation, dt = ku (k is the curvature), when we define 
the normal vector n = (ni, ri2, n 3 ) as 

n = n^i = -$ _ Vi$. (12) 

Thus, the normal vector n is given by the rotation of i (the unit vector along the x-axis). It 
is clear that t and n are orthogonal to each other. Now, 

dn = dniOi = -d($ - Vi$) = $-Vi$] = -k$"V 3 $ + r$"V 2 $, (13) 

which gives another Serret-Frenet equation dn = rb — nt, when we take the binormal vector 
b = {biMM) as 

b = b t a t = $-V 2 $. (14) 
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In this case, the binominal vector b is given by the rotation of j (the unit vector along the 
y-axis). Finally, similar calculation gives 

Ob = dbiOi = <9($~V 2 $) = -[^d®, $-V 2 $] = r$-Vi$, (15) 

which gives another Serret-Frenet equation, (9b = — rn. 
The time dependence of the r is similarly constructed as 

Br = -2i$- 1 a$$- 1 — $ + 2i$" 1 — (0$) = -2i$- x — — $ 

or or or 

= $ _1 (3r 2 o- 3 + 1kto 2 - -k 2 <j 3 + 9«<7i)$ 

= 3r 2 t + 2nrb - ]-n 2 t - dnn, (16) 

where t and b are given by Eqs. (10), (12) and (14). It is just the Eq. (3) and is already 
found by Lamb in [3]. This equation can be rewritten in terms of the curve variable r, 

Br = -d 3 r + (Jl<9 2 if - 3r 2 ) dr + 3rd 2 r x dr. (17) 

2.4 Equation for localized induction approximation with axial ve- 
locity 

Eq. (17) is the curve equation induced from the mKdV equation with the spectral parame- 
ter r. This equation can be rewritten to the equation for localized induction approximation 
generalized to take account of the axial-flow effect. Consider the following change of coordi- 
nates, 

z-+X-jT, z^T, r^-| (18) 
Above coordinates change transforms Eq. (17) into the following form, 

d T r = -8 x r - ^\d x r\ 2 d x r - cd\r x d x r, (19) 

where Ot = -§f,, dx = -§x- Eq. (19) is the equation introduced in [11] and it describes the 
motion of a thin vortex filament with axial velocity, embedded in a inviscid incompressible 
fluid. In [11], solutions of Eq. (19) were introduced including the N-soliton solution, a 
circular helix and a plane curve of Euler's elastica. But the Kida-type solution was not given 
explicitly and only equations describing the Kida-type curve was presented. Our results can 
be used to fill up this gap by using the transformation in Eq. (18). 
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2.5 Integrable structure of the curve equation 

We now use a modified form of the R-transformation [25, 20, 21] to obtain the Lax pair of 
the curve equation. It will be used to obtain a quasi-periodic solution by applying the finite 
integration method in the next section. First we define ^ = $~ 1( l> where $, $ are solutions 
of the linear equation (6) with the spectral parameter r, f, respectively . Define 

M = d^- 1 = + $- 1 <9<S<|- 1 <I> = $ _1 LW(t)$ - $ _1 LW(t)$ 

% 1 

= ^(r-T)^- 1 a 3 ^ = -(r-r)dr, 

2 2 V k / 



+ 4 (T " r)$ U« n 2 ) 



= % -(-2t 3 + 3fr 2 - r 3 )dr + \{t 2 - 2fr + r 2 )[dr, d 2 r] + Ut - f)Br. (20) 

Zi 4r _ 

Note that M 1 ]/ = <9^/, A^\l/ = '. Now, we introduce A = f — r, which becomes a new 
spectral parameter of the following linear equation for r, 

(9-M)* = 0, (<9-7V)* = 0, (21) 



where 



M = --A9r, 
2 



Q 1 

AT = - l -\ 3 dr - d\ 2 rdr + -X 2 [dr, d 2 r] - '-Xdr. (22) 

Zi Zi "rr _ 

Note that the compatibility condition of the associated linear equation, i.e., [d — M, d — 
N] = gives the following equation of motion for r, 

i-rd 2 r - ^[dr, d 3 r] - ^[dr, dr] = 0, (23) 

as well as the identity Tr(<9r) 2 = 2. Eq. (23) gives n and b components of dr in Eq. (2), 
which are U = —8k,V = 2kt in Eq. (3). But Eq. (23) does not give the t components 
of dr, i.e., W = 3r 2 — k, 2 /2 in Eq. (3). This fact is not surprising since W is known not 
to be related with the intrinsic form of the space curve, and is related with the internal 
parametrization of the curve [10]. In this respect, Eq. (23) just describes the intrinsic 
structure of the curve. In the next section, we will derive the quasi-periodic solution of the 
mKdV-type curve, starting from the associated linear equation in Eq. (21). In the course 
of derivation, we will additionally require the condition W = 3r 2 — k 2 /2, as well as the 
constancy of r, see section 3.4. 
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3 One-phase quasi-periodic solution 

3.1 Squared wavefunctions 

First we introduce notations, 

where = ctg = 0:3, /3+ = /3g = ^3,7+ = —7-, 73 = —73. The equation of motion 
in Eq. (23) is rewritten as 

Qirda 3 — ^73 = a + f3- — ot-(3 + , 

3irda + - = "3/3+ - a + f3 3 , (25) 

while the constraint Tr(<9r) 2 = 2 becomes a 2 + = 1. 

The finite-band integration method relies on the fact that the curve satisfying Eq. (23) 
guarantees the existence of the solution \l/ of Eq. (21) as a function of A. The following 
procedure assumes that the squared function of ^ is a polynomial function in A, and thus is 
called by the "modified squared wavefunctions method" [16]. Let the systems (21) have two 
basic solutions of column matrices, = (1^1,1^2) and ^2 = (<f>i, 02), which is used to build 
a squared wave function such that 

/ = -(i/2)(ip 1 (f)2 + ^20i), g = Mi, h = -4> 2 (j) 2 . (26) 

From the definition of /, g and h in Eq. (26) and using Eqs. (21), (22) and (24), we can 
obtain following equations, 

df = ~<^-9 + ^ a +h, dg = Xa+f - iXa 3 g, 
df = {\x'a + + \x 2 ra + + \x 2 1+ + \x^ + )h-{^ 

Bg = (-iX 3 a 3 -3iX 2 Ta 3 + ^X 2 -f 3 -iX[3 3 )g + (X 3 a + + 3X 2 Ta + + ^X 2 -f + + X[3 + )f. (27) 

Equations for h are similarly given. 

Using Eqs. (27), it can be explicitly checked that P(X) = f 2 — gh is independent of z 
and z and is only the function of A. The so-called 1-phase periodic solution is obtained by 
specially taking the form of 

P(X) = f 2 -gh = JJ (A - Aj), (28) 

i=l,4 
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where Aj are zeros of the polynomial which characterize the periodic solution of the "Bloch 
wave" problem and are called the main spectra. The zeros Aj have to consist of complex 
conjugate pairs Xj = Xrj + iXij, A J+ 2 = Xrj — iXjj, j = 1, 2, which is required for that the 
obtained curve r is real. Eqs. (27) and (28) are consistent with the following polynomial 
form of /, g and h in A, 

/ = a 3 X 2 + fiX + f , g = -iXa+(X- //), h = -iXa-(X - //), (29) 

where // are functions of z,z. 

3.2 Equations for \i, fa, fy, 7^ 

Using Eqs. (29) and (27), we can evaluate dg at A = /i to find 

( 9/x = -*(a 3 /i 2 + /iA' + /o)- (30) 
If we look at 0(A'),i = 2, 1, terms for df in Eq. (27), we find 

da 3 = --a + a^{fi - //*), = 9/ = 0. (31) 

At 0(A 2 ) terms for dg, we find 

da + = ia 3 a + fi + ia+fa. (32) 

Applying the same procedure to dg and df in Eq. (27), we can obtain that 

— % 
ifxa+d/i = (/i 3 a + + 3r/i 2 a + + -/i 2 7 + + /i/5+)(a 3 /i 2 + /i/i + /o) 

da + = ip 3 a + fi - ^7+/ + ip+fi + 3ira + f , 
= a 3 a + fi + a + f\ + -(7+^3 -73"+), 

1 1 

= -73«+At - f3 3 at+ + 3ra 3 a+/i + -7+/1 + «+/o + (3+a 3 + Sra+fa, (33) 



and 



— z — — 

da 3 = -(/3+a-ii* - P-a+fi), dfa = df = 0, 

= 2a + o;_(/i* — /i) — i(7 + o;_ — 7_a + ), 

= i(7 + a_ / u* — 7_a + /i) — 2(/3 + a_ — j3-Ct + ) + 6ro; + a_( / u* — //). (34) 
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3.3 n,fi,ji,(3i in terms of 

To prevent the "effectivization" problem, we start with Eqs. (28) and (29) to solve /i, f\ and 
f in terms of a 3 . First, we introduce constants of motion Sj, i = 1,4 which are defined 

P(A) = f - gh = A 4 - siA 3 + s 2 A 2 - s 3 A + s 4 . (35) 

By inserting /, (7 and /i in Eq. (29) into Eq. (35), we can obtain 

si = J^Ai = -2/ia 3 + a+a;_(/j + //), 
S2 = 51 = 2 /o«3 + fi + a+a-iifj,*, 

i<j 

s s = AiAjAfc = -2/1/0, s 4 = n^* = /o- ( 36 ) 

i<j<k 

Now solving for / , / 1? /x and /i* in Eq. (36), we obtain 

/1 = ttt^, /o = -v^, V, V* = 7771 — 2\ (Vs^az + Si-s 4 =F JhR] , (37) 

2V S 4 A 1 ~ «3J S 4 V / 

where 

R = (x + A 2 A 3 + AiA 4 )(x + A1A3 + A 2 A 4 )(:r + A X A 2 + A 3 A 4 ), 

x = 2 v / s 4 a;3. (38) 

Note that f i: i = 0, 1 are constants, which are consistent with Eqs. (31) and (34). To obtain 
this solution, we use Mathematica for symbolic manipulation. In addition, Mathematica was 
used to check various formulae appeared in this paper. For example, it was used to check 
that Eqs. (55) and (64) satisfies the equation of motion (23). 

The definition of 7^,2 = 1,3 in Eq. (24) together with Eqs. (31), (32) and (36) give 

73 = -ia + a_([j, + n* + 2/ 1 /a 3 )a 3 = -i(2fi/a 3 + si)a 3 = ta 3 , 

7 + = {2i(/i + 03/1) — ia + a_(fi + fi*)}a + = 2(2/1 — Si)a + = ya + , (39) 

and 7_ = —y*a^. Note that these results are consistent with the third equations of Eqs. 
(33) and (34). 

Now, we require that 0:3 becomes a traveling wave in the 1-phase solution, i.e., a function 
of w = z + vz, where v is a constant, requiring da 3 = vda 3 . The definition of a 3 , (5 3 in Eq. 
(24) gives da 3 = d(3 3 , which in turn becomes 

(3 3 = va 3 + c(z). (40) 
11 



To fix c(z), we need following steps. The fourth equations of Eqs. (33) and (34) together 
with the results in Eqs. (39), (40) give 

P+ = {fi/al + fif 1 /a 3 -fo/a 3 + v + c(z)/a 3 --(f 1 /a 3 + fi)t-3r(f 1 /a 3 + fi)}a + = ua + , (41) 

and P_ = u*a_. Now, the first equation of Eq. (34) becomes 

_ i 

da 3 = -a + a-(un* — u*fi). (42) 
Now, the relation da 3 = vda 3 with Eqs. (31) and (41) gives c(z) = c is a constant and 

c = f + fi8 1 /2 + 3rf 1 . (43) 
Finally, inserting all these results into the first equation of Eq. (33), we can obtain 

dfj, = -i(fi 2 + 3r/i + -fiy + u)(a 3 fi 2 + frfi + / ) = vd/i, (44) 
which shows /i is also a function of w = z + v z only. 

3.4 Imposing the conditions W = 3r 2 — k 2 /2 and r = constant 

To find the value of v, we need to impose the condition W = 3r 2 — /t 2 /2, i.e., Trdrdr = 
6r 2 — k 2 . First, we calculate the curvature k. Using Eqs. (11), (24), (31), (32), (36), we find 

2k 2 = Tr(d 2 r) 2 = -Lsj + 2s 2 - Aa 3 f . (45) 

Similarly, we can obtain 

Trdrdr = 2v + 2a 3 f - sj/2, (46) 
which together with Eq. (45) gives 

v = -sI-^s 2 + 3t 2 + -t Si . (47) 

Now, we calculate the torsion and impose the condition that the torsion is a constant. 
We first find {t, n, b} vectors as following, 



ha; = dr 



a 3 a + 
a_ —a 3 



n . a . = 1^. (j. — — ( a - a +(^* ~ f 1 )/ 2 a+{a^ + h) 

k 1 1 k V -a-(a 3 /i* + fi) -a-a + (n* - n)/2 
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Inserting the result in Eq. (48) into the equation <% 3 = — rn 3 gives (use Eqs. (45) and (31)) 

(y + ^ 2 + /o(y«3 + /i) = 0, (49) 

which can be satisfied with a constant r by choosing (use Eq. (45)) 

si = -4r, s 3 = 2r(4r 2 - s 2 ). (50) 

We check that the curvature k given by Eq. (45) becomes the solution of the mKdV equation 
when si, s 3 satisfies the conditions in Eq. (50). Using Eq. (36), the condition in Eq. (50) can 
be expressed in terms of X R ±, X R2 , An, A/ 2 where Aj = X R j + iXij, Xj +2 = X R j — iXjj, j = 1, 2. 



type-A : X R1 = X R2 = — r, An and A/ 2 are arbitrary, 
type-B : X R1 = -X R2 - 2r, A n = X I2 . 



Thus there exist two types of curves of the mKdV equation, which require different choice 
of the main spectra Aj. Each type of main spectra has two independent parameters for curves 
of torsion r. They are An and A/ 2 for the Type-A, while they are X R \ and An for type-B. 
In addition, we require the condition on v in Eq. (47). 

3.5 Derivation of r-j 

Using all these results, we can obtain the equation for a 3 , 

dot3 i 



4. 2vv A (51) 

This equation can be integrated in terms of Weierstrass' V(u,g 2 ,g 3 ) function. As far as 
Weierstrass elliptic functions are involved, we employ terminology and notation of [26] with- 
out further explanations. Explicitly 

a 3 = ^= (-j- 4V(w + w 3 , Q2, <? 3 )) , (52) 

where w 3 is an integration constant and 

1 2 1 

92 = —S 2 - -S1S3 + s 4 , 

1 3 1 1 2 1 1 2 

93 = ^Q-S 2 - -^S lS2 S 3 + — S 3 - -S 2 S 4 + T^ S 1 S4 - ( 53 ) 
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The integration constant w 3 is determined by the initial condition, which we shall choose as 
follows; V(w 3 ) = e 3 at w = 0, where e 3 is the smallest root of the equation 4z 3 — g 2 z — g 3 = 0. 
This condition guarantees that | ck 3 j < 1. Other two roots are denoted by e\ and e 2 with 
ei > e 2 > e 3 . w 3 as well as w\ are called the half period of the V function. They satisfy 
V(wi) = ei, V(w 2 ) = e 2 , e\ + e 2 + e 3 = w 1 + w 2 + w 3 = 0. Especially, 

ei = ^(AiA 3 + A 2 A 4 ) - ||, e 2 = ^(AiA 4 + A 2 A 3 ) - ||, e 3 = ^(AiA 2 + A 3 A 4 ) - ||. (54) 

Now using a 3 = <9r 3 and /3 3 = (9r 3 in Eq. (40), we can obtain (use Eq. (43)) 

s 3 si - 4s 4 + 6rs 3 _ 



r 3 = y a 3 d-u; + 



2 s 2 1 

— {((w + w 3 , g 2 , g 3 ) - -|w + -(s 3 si - 4s 4 + 6ts 3 )z - r] 3 }, (55) 

S 4 12 o 



where ({u,g 2 ,g 3 ) is the Weierstrass' zeta function and the integration constant r] 3 is taken 
to be ((w 3 , g 2 , g 3 ). (We will use the notation rji = ((wi, g 2 ,g 3 ),i = 1,3 in the following.) 

3.6 Derivation of r + 

Now we try to obtain the solution r + = r\ — ir 2 of the mKdV-type curve. The procedure is 
similar as in [20, 21], and we will only describe some important steps of derivations. More 
details including the identities of Weierstrass' functions needed in the derivation can be 
found in [20, 21]. Using Eqs. (39), (40), (41), we can express the second equation of Eq. 
(33) as 

da + = vda + + i{s\/2 + 3r)f a + , (56) 
which, with the help of Eq. (32), gives 

a + = exp{— iy / s 4 (si/2 + Sr)z} a + (w), (57) 

where a + satisfies the following differential equation, 

^± = ia + (f\ + a 3 /i). (58) 

Using Eqs. (37), (51), (52), we can integrate Eq. (58) to obtain (see more details in [20, 21]) 

o * <x(w + w 3 + K^ajw + w 3 + K 2 ) 

a+ = 2—= r exp{-C(fi;i)w -Q(k 2 )w + 6\, (59) 

a/s 4 a{Ki)a{K 2 )a 2 {w + w 3 ) 
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where two constants k±,k 2 are defined by the following relations (Please distinguish the 
constants K\ and k 2 from the curvature k.), 



V(k 2 ) = -s 2 /l2 ± (60) 
V f ( Kl ) 1 V(K 2 ) = -^^\ =~VR\ =~(siVi4=FSs). (61) 

2 (IW a 3 =Tl 4 a 3 =Tl 4 

The integration constant 5 can be fixed by requiring a 3 + |«+| 2 = 1. Explicit calculation of 
this requirement, especially at u> = 0, gives 5 = —t] 3 (ki + k^)- 
We now derive r + . Using dr + = a + in Eq. (24), we obtain 

r+ = J a + dz + M(z) = exp{ -iy/sl(si/2 + 3r)^} J a+dw + M(z), (62) 

where M(z) is a function to be determined. To evaluate the integration explicitly as well as 
to determine M(z), we substitute Eq. (62) into the following relation, 

dr + — j3 + — ua + = (v — si/i/2 — 3T/i)a + . (63) 

It then gives M(z) = and (see more details in [20, 21]) 
2i 

r+ = —{((w + W 3 + K 1 )-((w + W 3 + K, 2 )-(( y K, 1 )+((K 2 )} 
S4 

a(w + w 3 + Ki)a(w + W 3 + K 2 ) 
a(K 1 )a(K 2 )a 2 (w + w 3 ) 
x exp{-C(/ti)-u; - ((k 2 )w - tj 3 (k,i + k 2 ) - i y /s2(s 1 /2 + 3t)z}. (64) 

Eqs. (55), (64) are the main results of the present paper. The space curves of the mKdV 
equation are described by taking two types of the main spectra \,i = 1,4 shown in section 
3.4. They are curves which have a constant torsion r and their curvature k is given by Eqs. 
(45) and (52). 

4 Special Cases 

In this section, we study some special cases of the obtained solution by taking specific values 
on the main spectra A,. The Weierstrass functions reduce to simple forms (for example, 
sinusoidal or elliptic functions) on these special limits. These cases contain the straight 
line, the Kelvin wave and the 1-solitonic curve, as well as plane curves in terms of Jacobi 
functions, and the closed ring. 
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4.1 The straight line 

The simplest solution that corresponds to the straight line is obtained by taking An = A/2 = 
in the type-A curve (or \ R i = — r, An = in the type-B curve), i.e., Ai = A 2 = A 3 = 
A 4 = — t (real). In this case, g 2 = g 3 = e 3 = 0. The Weierstrass functions in this limit 
are given by V(u,0,0) = l/u 2 , ((u,0,0) = 1/u, a(u, 0, 0) = u. And w 3 — k± — oo. Using 
these relations, we can easily obtain r 3 = —z + 3r 2 z,r + — r 1 — ir 2 = 0. The rotational and 
translational symmetry of the vortex equation can give more general configuration of the 
straight line [20]. 

4.2 The Kelvin wave 

The next example, known as the Kelvin wave in fluid mechanics [27], corresponds to taking 
A/i = in the type-A curve, i.e., Ai = A3 = — r, \ 2 = — r + ry, A4 = — r — ry. In this case, 
92 = j2l 4 ,93 = 2l67 6 - As A = (?| — 277! = 0' tli e Weierstrass functions are given in simple 
forms, 

V(u) = - 7 2 /12 + 7 2 csc 2 (7m/2)/4, 
£( u ) = 7 2 M/i2 + 7cot(7w/2)/2, 

a(u) = 2exp(7 2 M 2 / 24 )sin(7M/2)/7. (65) 



And e 2 = e 3 = — 7 2 /12, w 3 = ioo, 

d\r\( rvKTr, l0\ = 

y^T 2 ± 2T^fT T T^i 



sin(7/ci/2),sin(7/c 2 /2) = ; , (66) 



cot(7«i/2), cot(7/t 2 /2) = -- ( ^/r 2 + 7 2 ± r ) . (67) 



Now, a straightforward calculation gives 

1 



r 3 = ^— ^ {^ + 3(r 7 2 /2 + r 3 )^}, 

J_ 

T 2 + -f' Z 



r+ = i - 2 exp{i^r 2 + 7 2 (^ + t 2 z - Yz/2)}. (68) 



This curve is a helix (r + = r 1 — ir 2 ). When r = 0, it describes a circular ring having a 
constant radius and rotating on the x|/-plane. The parameter 7 is the curvature n of the 
curve, i.e., k — 7, while r is the torsion of the curve. The three orthogonal vectors are given 
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by 

t = . =(—k cos Q, k sin Q, r), 
n = (sin 6, cos 6,0), 

b = =(r cos 9, — r sin 6, k), (69) 

yr + k 

where 6 = \Jt 2 + /t 2 (,2 + r 2 z — k 2 z/2). 

4.3 The Hasimoto 1-solitonic curve 

To obtain the curves of 1-soliton from our periodic solution, we take An = A/2 in the type-A 
curve (or A#i = — r in the type-B curve), i.e., Ai = A2 = — r + 27, A3 = A4 = —r — ry. 
In this case # 2 = fT 4 ,^ = -J^ 6 and w 3 = in/(2~f), e 3 = P(tu 3 ) = -27 2 /3, Vs = CM = 
— ?77r/6, <t{wz) = iexp(7r 2 /24)/7. The Weierstrass functions are given by 

V{w) = -27 2 /3 + 7 2 coth 2 (7w), 
C(w) = — 7 2 w/3 + 7 coth(7w), 

<t(w) = exp(— 7 2 w 2 /6) sinh(7w)/7. (70) 

And = 00, sinh7/t 2 = —ij/^r 2 + 7 2 , coth7/t 2 = rr/7. It then gives 

27 

r 3 = —z — 3r 2 z + — tanh(7^ - 7 3 z), 



T 2 + Y 



2«7 

r + = -sech(72 — 7 3 z) exp(— irz — ir 3 z). (71) 



The curvature of the curve is k — 27sech(7^ — 7 3 z). It satisfies the mKdV equation in 
Eq. (7). The three orthogonal vectors can be easily constructed. For example, 
27sechi~2 

t = — — (7sin6tanhf2 — tcos6, — 7cos6tanhf2 — rsinG, 

T Z + Y 

r 2 +-y 2 

7 sechfi cosh Q) , (72) 

where Q = 72: — 7 3 z, O = rz + r 3 z. In Figure 1, we show an example of the 1-solitonic curve 
with parameters 7 = 1/3 and (a) r = 1, (b) r = 3. 

4.4 The rigid vortex filament 

When we take A/ 2 = A^ 2 + t in the type-B curve, i.e., Ai = A 3 = — a — 2r + i(a + r), A 2 = 
A4 = a + i(a + r) for arbitrary a, then v — (use Eq. (47) ) and w = (z — vz) — > z. Thus 
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Figure 1: Typical 1-soliton curves drawn for —10 < z < 10, (a) with r = 1, 
(b) with r = 3. The parameter of the curves is 7 = 1/3, drawn at z — 0. 

it describes a rigid filament of fixed shape, which rotates around the z-axis with constant 



angular velocity (=— rw4(a + r) 4 + r 4 ) and moves along the z-axis with constant velocity 



(=— (4(a + r) 4 + 3r 4 )/y 4(a + r) 4 + r 4 ). This configuration could be easily identified in 
experiments. 

4.5 Plane curves in terms of Jacobi functions 

Another reduced form of our solution is obtained by choosing K\ = w 1 , k 2 = w 3 . Then 
V(ki) = V(uo\) = ei,V(K 2 ) = V(u 3 ) = e 3 , and Eqs. (60) and (54) fix the main spectra as 
Ai = A3 = pexp(i6), A 2 = A4 = — pexp(— iO). Note that it corresponds to taking r = in 
the type-B curve, i.e., \r\ = pcosO, An = psin6> . The Weierstrass a function is reduced to 
the Jacobi's elliptic function under this limit as following, 
a(w + w 3 + Ki)a(w + w 3 + k 2 ) 



a 2 (w + w 3 ) 



exp (^—{w 2 + w 2 — 2w 2 w 3 — 2w 2 w) ) 
V 2vj\ J 



xq 3/4 exp ( - *2^r w ) \j ^ sn ^ ei ~ e 3 w )^ n W e i ~ e 3^). (73) 



Using the relation r)iw 2 = — in/ 2 + i] 2 wi (Legendre's relation) and k = \Je 2 — e 3 je\ — e 3 = 
sin6*, k' = cos 9, q = exp(inw 3 /wi), we can obtain 

r + = 2isin9 cn(pw)/ p, 

r 3 = J (-l + 2dn 2 (pw)^)dw-p 2 z = -w + 2E(pw\k)/p-p 2 z, (74) 
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Figure 2: (a) Open plane curve for —10 < z < 10 and z — 0. (b) Closed 
plane curve showing the famous figure-8 shape. 

where E is the incomplete elliptic integral of the second kind. The curve lies on the yz- 
plane, and has the curvature k = 2kpcn{z — p 2 (2k 2 — l)z} and torsion r = 0. Figure 2 
show examples of the plane curve. Fig. 2(a) shows an open plane curve which is obtained 
by taking Ai = A* = (3 + 3i)/v% A 2 = A* = (-3 + 3i)/y/2,r = 0. Fig. 2(b) shows the 
famous figure-8 shape (closed plane curve) in [4, 10], which is obtained by taking special 
value 6 — 1.14. The condition for closed curves will be discussed in the following section. 

4.6 Closed curves 

The characteristics of closed curves can be obtained from the following quasi-periodicity 
properties of Weierstrass' elliptic functions; 

t{u + 2wi) = ((u) + 2r ]i , 

a(u + 2wi) = — exp{2i]i(u + Wi)}a(u). (75) 

In our problem, W\ is real, while is pure imaginary. Thus the physical characteristics of 
closed curves are described by the period Aw = 2w\. After this period 2w\, r + obtains an 
additional factor 

exp{2?7i(/ti + k 2 ) - 2((k 1 )w 1 - 2((k 2 )w 1 }. (76) 
Thus a necessary condition for a closed curve , i.e., r + (z = 2mwi) = r + (z = 0) , is 

771 (ki + k 2 ) - iui{C(«i) + ((^2)} = i— 7T, (77) 
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Figure 3: A typical closed ring plotted for — 10 < z < 10 with parameters 
\ 1 = A* = -2.0653 + 1.17H, A 2 = A* = 0.0653 + 1.17H, r = 1 (a) at z = 0, 
(b) at z — 5, (c) at z — 10.09 (quasi-period time). 

where m, n are arbitrary integers. Another condition comes from the closedness of the z- 
coordinate, i.e., r 3 (z = 2mwi) = r 3 (z = 0), which is 



52 1 n 

-Wi + —^zT)i = 0, 



(78) 



i.e., s 2 wi = 12771. 

The quasi-periodic property of the closed curves along the time is following. During a 
time period Az = 2w 1 /v, the ring returns to its original shape, but rotates around the z-axis 
by an angle 

2-7r- v /i7( Sl + 6rA (79) 
m v 

and moves a distance 

S3S1 - 4s 4 + 6rs 3 wi , OAN 

(80) 



2 v /s 4 " V 
along the z-axis. 

The two conditions, Eqs. (77) and (78), can be solved numerically using a software 
package like Mathematica. One possible solution is Ai = A3 = —2.0653 + 1.17U, \ 2 = \* 4 = 
0.0653 + 1.171i,r = 1, which belongs to type-B curves. It satisfies the conditions (77) and 
(78) with n = 2,m = 5 (w\ = — 1.198, v = —0.2374). Figure 3 shows the motion of a closed 
curve at (a) z — 0, (b) z — 5, (c) z = 10.09 (—Az, quasi-period time), respectively. During 
a time period Az, the ring rotates -25.6 radian and advances -37.07 along the z-axis. The 
z period is 10u>i = —11.98. We can't find the type-A curve that satisfies the closedness 
conditions in Eqs. (77) and (78). 

In the case of figure-8 shape (closed plane curve), the condition K\ — w±, k 2 — w 3 satisfies 
Eq. (77) with n — l,m — 2 (Legendre's relation.). The second condition (78) becomes 
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Figure 4: Typical open curves having r = 1, drawn for — 10 < z < 10, z = 0, 
(a) type-A curve with Ai = A3 = — 1 + i, A 2 = A4 = — 1 + 2i, and (b) type-B 
curve with Ai = A3 = -3 + 2i, A 2 = A4 = 1 + 2i. 

G((wi) = —p 2 w 1 , which is satisfied by choosing Ai = A3 = p(l + 2.177966i), A 2 = A4 = 
—p(l — 2.177966i), or 9 = 1.14 in the notation of section 4.5. Fig. 2(b) shows the figure-8 
shape curve. 

4.7 General vortex filament 

Most generally, the solution in Eqs. (55) and (64) describes an open space curve with a 
constant torsion r. Figure 4(a) shows an example (type-A curve), which we obtain using 
Ai = A3 = — 1 + i, A 2 = A| = — 1 + 2i, t — 1. Figure 4(b) shows another example (type-B 
curve), which we obtain using Ai = A3 = — 3 + 2i, A 2 = A4 = l + 2i, r = 1. These are examples 
of the most general configurations described by the 1-phase quasi-periodic solution. 

5 Discussions 

In this paper, we present the calculation of the one-phase quasi-periodic solution of space 
curves described by the mKdV equation. The solution is explicitly given in terms of Weier- 
strass' elliptic functions. The main difference of the present calculation compared to that of 
previous works in [19, 20, 21] is the appearance of constraints on the main spectra resulting 
from imposing two conditions. It then results in two types of curves, type-A and type-B 
curves. The first condition, W = 3r 2 — k 2 /2 in Section 2.4, is related with the internal 
parametrization of the curve. It does not determine the extrinsic form of the curve [10]. 
This is consistent with the fact that the curve equation in Equation (23) only fixes the form 
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of U and V and does not fix W in Equation. (2). As is shown in Section 3.4, the condition 
on W fixes the velocity v of the motion of rigid curve. The second condition, r = constant, 
is also unique in the present calculation, and is not required in the previous calculations in 
[19, 20, 21]. It was checked that the constraint on the main spectra is required for that k in 
Equation (45) satisfies the mKdV equation in Equation (7). 

We study various configurations resulting from the degenerate limit of the quasi-periodic 
solution. They contain the Kelvin wave belonging to the type-A curve, the rigid vortex 
filament, the Jacobian plane curve, and closed curves belong to the type-B curve, and the 
straight line and the one-solitonic curve belonging both to the type-A and type-B curve 
simultaneously. 

The simple case of zero-torsion was studied in various contexts such as the elastica prob- 
lem and the computer vision problem. In this case, the fact that the curve lies on a plane 
was the main keystone in the derivation of the curve configurations. In the case of non-zero 
torsion, the convenient formalism of differential geometry seems difficult to realize except 
in some especially simple cases. The present derivation heavily depends on the integrability 
structure of the mKdV equation. The equations for the one-phase solution in Section 3.2 are 
by no means simpler that the equations for the stationary curve configuration given in [11] 
(Equations (3.32), (3.36), (3.38), and (3.39)). But there exists a standard procedure in solv- 
ing these equations, shown in Section 3.2. Moreover, this procedure offers a generalization 
which can be used in calculating the A-phase periodic solutions. 

Our formalism might be adapted to the curve problem of the non-zero and non-constant 
torsion problem. It is known that the spectral parameter of the integrable equation can be 
extended to be z, z-dependent without destroying the integrability [28, 29]. Thus it could 
result in space curves of non-constant torsion. Explicit construction of curves with non- 
constant torsion is rare and difficult. Thus it remains as an interesting work for the future. 
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